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In this paper, linear first order expansion of deceleration parameter q{z) — qo + qi(l — a) (Mi), 
constant jerk j = jo (M2) and third order expansion of luminosity distance (M3) are confronted 
with cosmic observations: SCP 307 SN la, BAO and observational Hubble data (OHD). Likelihood 
is implemented to find the best fit model parameters. All these models give the same prediction of 
I the evolution of the universe which is undergoing accelerated expansion currently and experiences a 

transition from decelerated expansion to accelerated expansion. But, the transition redshift depends 
on the concrete parameterized form of the model assumed. Mi and AI2 give value of transition 
^2 ' redshift about zt ~ 0.6. M3 gives a larger one, say zt ~ 1. The /d-of implies almost the same 

, goodness of the models. But, for its badness of evolution of deceleration parameter at high redshift 

' z > 1, Ms can not be reliable. Mi and M2 are compatible with ACDM model at the 2a and 

00 ' la confidence levels respectively. M3 is not compatible with ACDM model at 2a confidence level. 

(N ■ From Ml and M2 models, one can conclude that the cosmic data favor a cosmological model having 

I ^- jo<-l. 
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■ I. INTRODUCTION 

1^1 ' The expansion of the universe or kinematics of the universe is described by the expansion rate H = d/a, a 
dimensionless parameter q — — aa/a^, named deceleration parameter and jerk parameter j — — da^/[aaP), where a 
T-H ■ is the scale factor in Friedmann-Roberson- Walker (FRW) metric. The front negative sign of q is added to obtain a 
^ ' positive parameter when one considers a decelerated expansion universe dominated by matter fields with attractive 
force, for example dark matter dominated universe with q = 1/2. However, the observations of Type la Supernovae 
(SN la) from two teams [HQ imply that our universe is undergoing an accelerated expansion at present. Whereafter, 
this result is confirmed by the observations from WMAP and Large Scale Structure survey @. Then at 

present, the deceleration parameter would be a negative number. In fact, to explain the accelerated expansion of the 
universe, a fiood of various cosmological models have been explored, please see [T/] for recent reviews. 
^ ^ In the literatures, two approaches have been taken into accounts. One is that an extra energy component, dubbed 

[ dark energy, is introduced, e.g. cosmological constant, quintessence [1, [IS [HI j phantom and quintom [l^, etc. 
The other is that the accelerated expansion of the universe is due to the modification of the gravity theory at large 
scale, e.g. modified gravity theory, Brans-Dicke theory and higher dimensional theory, etc. These correspond to the 
dynamics of the universe. In general, different dark energy models would predict different expansion histories of the 
?H ■ universe, for example quintessence dark energy model with various scaler potentials will have different values of 
the transition redshift zt from decelerated expansion to accelerated expansion. These different models would have 
different dynamics. However, a kinematic approach will be held regardless of the underlying cosmic dynamics [l5l.[l^. 
In particular, the jerk parameter can provide us the simplest approach to search for departures from the cosmic 
concordance model, for its constant value j = —1 for cosmological constant. This approach is called cosmography 
[13, [Hi, cosmokinetics JJ.], or Friedmannless cosmology [20, 21]. 

Some explore the accelerated expansion of the universe by using different parameterized forms of deceleration 
parameters q{z) [2^[23l.[2^ in the so-called model independent way, for examples constant model q = constant, linear 
model with variable z {q{z) = qo + qiz) and linear model with variable a {q{a) = go + ~ o,)), etc. Recently, the 
authors of [l6l . [2l| investigated constraints on some kinematic models by employing a Bayesian marginal likelihood 
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analysis based on the cosmic observations from supernovae type la with extent including the jerk parameter j: the 
third order contribution in the expansion for kinematic luminosity distance in terms of the redshift z. 

As known, the degeneracies between the parameters will be broken when more different observational data added 
in constraint to models. In this work, we will go further by including BAO and the observational Hubble data in 
examining the case for the jerk parameter (constant and variable). 

II. KINEMATIC MODELS 

The metric of a flat FRW cosmological model containing dark matter and dark energy is written as 



(1) 



where a{t) is the scale factor, which describes the whole history of the universe evolution, and has the relations with 
redshift z in terms of a = (1 + z)~^ (ao = 1 is normalized). The Hubble parameter 



a 



and deceleration parameter 



1 a _ 1 , 



1. 



(2) 



(3) 



are defined as the rate of expansion and accelerated expansion. By using the relation oo/a = 1 + z and the relations 
between H and i.e., one can rewrite Eq. Q in its integration form 



H{z) ~ Ho exp 



[l + q{u)]d\n{l+u) 



(4) 



Similarly, the jerk parameter is defined as 



J 



1 d 



1, ,2[Hiz)^]' 



(1 + ^) 



[H{z 



1 



H{zY 

Easily, one can find that the deceleration parameter and jerk parameter have the relations 



(5) 
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q + 2q' + {l + z)^ 
dz 



(6) 



which will be used when the parameterized forms of q{z) are giong to be tested. By using these definitions Eq. 
Eq. ([3]) and Eq. ([5]), one can describe the recent cosmic expansion with their current values 



lit) = 1 + Hoit - to) - \qoHl{t - to)^ - ^JoH^it - to)3 + 0[it - to)^] 



from which the luminosity distance can be expanded as [2 



c 



^ + ^(1 - go)2^ - i(l -qo- iql - jo)z^ 



0{z^ 



(7) 



(8) 



Also, one can find the relations of Hubble parameter H{z) and luminosity distance dL{z) [2f 

H-\z) = -il + z)^ 
dz 

= ^Jil + z)-'dUz)] 



(9) 



which will be useful when the observational Hubble data are used as cosmic observation constraint. 

The basic aim in this paper is to examine some simple kinematic models for the cosmic expansion based on specific 
parameterizations for q(z) in Eq. ([3]) (variable jerk parameter), a constant jerk parameter and their comparison with 
the expansion ([7]). 
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The first model, Mi, is given by linear expansion of the scale factor a, q{a) = go + ~ a) in terms of a, which can 
be rewritten in the terms redshift z, q{z) = qQ + qiz/{l + z). The second model, M2, is a constant jerk parametrization 
model, j{z) — jo, for detecting the departure from the flat ACDM scenario, for which j(z) = jo = —1. Model M3 is 
the expansion ([7]), which has as free parameters go and jo- In the Appendix [XI see also [16j, one can find the basic 
analytical expressions between the Hubble parameter H{z), deceleration parameter q{z) and jerk parameter j{z). One 
can take this paper as a generalization and complement to pj6| where only the 307 SN la data points are used. Here, 
the BAO and OHD datastes are also included as useful cosmic constraints. One would notice that all of them do not 
include Q.m explicitly. 



III. COSMIC OBSERVATION DATA SETS AND STATISTICAL RESULTS 

A. SN la 

We constrain the parameters with the Supernovae Cosmology Project (SCP) Union sample including 307 SN la 
[2^ . which distributed over the redshift interval 0.015 < z < 1.551. Constraints from SN la can be obtained by fitting 
the distance modulus 



^ith{z) = 51ogio(L'L(z)) + Ato, 
where, Di^{z) is the Hubble free luminosity distance Hodi^{z)/c and 

dz' 

ddz) = cil + z)J^ ^ 
fio = 42.38- 51ogio/i, 



(10) 

(11) 
(12) 



where Hq is the Hubble constant which is denoted in a re-normalized quantity h defined as Hq — lOO/i km s ^Mpc ^ . 



The observed distance moduli ^obs{zi) of SN la at Zi is 

fJ-obs{zi) = mobsizi) - M, 



(13) 



where M is their absolute magnitudes. 

For SN la dataset, the best fit values of parameters in a model can be determined by the likelihood analysis is based 
on the calculation of 



X^{Ps,mo) = ^ 



[^J■obs{z^) - llth{Ps,Zi)y 



SNIa 



E 

SNIa 



[5logio{DL{ps,Zi)) -mobs{zi) +moY 



(14) 



where mo = /xg + M is a nuisance parameter (containing the absolute magnitude and Ho) that we analytically 
marginahze over [27| . 



/+00 
exp 
-00 



to obtain 



where 



£2 



X = A - — + — 



]pC^{Vs,mo) 



C 



dmo , 



C 



27r 



[5 \ogioiDL{ps,z,)) -■niobsizi)]'^ 



(15) 



(16) 



(17) 



SNIa 



B = 



5 logio (Z?!, {ps ,Zi) - ruobs {zi ) 



SNIa 



(18) 
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^=E^- (19) 

SNIa ' 

The Eq. (I14p has a minimum at the nuisance parameter value toq — B/C. Sometimes, the expression 

xlNiaiPs,B/C)=A-{ByC) (20) 

is used instead of Eq. ^TE\i to perform the hkehhood analysis. They are equivalent, when the prior for niQ is flat, as 
is implied in (jl5[) . and the errors ct^ are model independent, what also is the case here. Obviously, from the value 
mo = B/C, one can obtain the best-fit value of h when M is known. 

To determine the best fit parameters for each model, we minimize x^iPs, B/C) which is equivalent to maximizing 
the likelihood 

Cips) oc e~x'(^'-^/^)/^ (21) 
B. BAO 

The BAO are detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples, and measure 
the distance-redshift relation at z 0.2. BAO in the clustering of the SDSS luminous red galaxies measure the 
distance-redshift relation at 2; = 0.35. The observed scale of the BAO calculated from these samples and from the 
combined sample are jointly analyzed using estimates of the correlated errors, to constrain the form of the distance 
measure Dyiz) [H, [ll [s^l 



Dviz) 



1/3 



H{z) 

where Da{z) is the proper (not comoving) angular diameter distance which has the following relation with (ii(z 



(22) 



D.i,) = (23) 

Matching the BAO to have the same measured scale at all redshifts then gives [30| 

£iy(0.35)/Dv(0.2) = 1.812 ± 0.060. (24) 

Then, the Xbao(Ps) is given as 

2 , , [Dv{0.?,b)/Dv{0.2) - 1.812]^ 
Xbao{Ps) = • (25) 

C. OHD 

The observational Hubble data are based on differential ages of the galaxies Isil. In [3^ , Jimenez et al. obtained an 
independent estimate for the Hubble parameter using the method developed in |3l|, and used it to constrain the EOS 
of dark energy. The Hubble parameter depending on the differential ages as a function of redshift z can be written in 
the form of 

H{z) = --^^ (26) 
1 + z at 

So, once dz/dt is known, H{z) is obtained directly [11]. By using the differential ages of passively-evolving galaxies 
from the Gemini Deep Deep Survey (GDDS) [13] and archival data [ssl . [s^JstI. [ssl [39l . l40j . Simon et al. obtained 
H{z) in the range of < z < 1.8 [Ij]. The observational Hubble data from^sjare list in Table HI 

The best fit values of the model parameters from observational Hubble data [ssj are determined by minimizing 

2 / N _ V" [Hth{Ps]z,) - Hobs{zi)f 
XHub{Ps)-y^ -T^^ , (27) 

where Ps denotes the parameters contained in the model, Htu is the predicted value for the Hubble parameter, -ffobs 
is the observed value, a{zi) is the standard deviation measurement uncertainty, and the summation is over the 9 
observational Hubble data points at redshifts Zi. 
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z 


0.09 


0.17 


0.27 0.40 


0.88 


1.30 


1.43 


1.53 


1.75 


H{z) (km s-^Mpc"^) 


69 


83 


70 87 


117 


168 


177 


140 


202 


la uncertainty 


±12 


±8.3 


±14 ±17.4 


±23.4 


±13.4 


±14.2 


±14 


±40.4 



TABLE I: The observational H{z) data [3i,|3l|. 



D. Statistical Results 

For Gaussian distributed measurements, the likehhood function L oc e"^ where is 

= XSNIa + xIaO + XHub ' (28) 

where XsNia^ x'bao ^^'^ X^Hub ^""^ ^^^^ described in Eq. (PO)) . Eq. (pS)) and Eq. P7)) respectively. It is clear 
that only the kinematic variables, say /lo, 90 and jo, are contained in all x^s equations where ilm does not appear 
explicitly. Then, the results may not depend on the dynamic variables, say and gravitation theory. Here more 
datasets are included to constrain the model parameters than that in Ref. jl6| where the SCP 307 is used alone. 
After the calculation as described above, the result is listed in Table ITTl 



Models 


Xmin 


50(10-) 


JO (la) 






xVdof 


Ml 


326.599 


-0.715±roS 


Z.iyD_o 241 


716+0 053 
J^O-0.053 


0.609±«;J^° 


1.043 


Mz 


326.442 


-0.658±S;;«^ 


1 009+0.219 

1.00Z_g 225 


7nq+0 053 

U. (Uy_o 053 


0.592tH^6! 


1.040 


Ms 


330.733 


-0.46llH^3^ 


U-J^*' -0.094 


716+0 053 


0.946l«;J^? 


1.053 



TABLE IL The calculation results of the models. The dof denotes the degrees of freedom of the models. 

It can be easily seen that Mi and M2 almost give the same result and have the same goodness in the viewpoint of 
/dof . There the jerk parameter jo is different: Mi has larger absolute value than that of M2. However, M3 has 
smaller absolute value of go and ja than that of Mi and M2 and predict larger transition redshift from decelerated 
expansion to accelerated expansion. The corresponding evolution curves of deceleration parameter q{z) with respect 
to redshift z are plotted in Fig. [T] From the Fig. [l] one can see that Mi and M2 almost give the same evolution 
history of the universe. The difference between them is that the q{z) of M2 is depressed at high redshift than that 
of Ml. However, M^ almost gives the wrong evolution history of the universe at high redshift for q approximates to 
0.5 when the epoch of matter dominated. The reason is simple that the model M3, expansion of scale factor at low 
redshift, is not reliable at high redshift. But, we must notice that M3 is the worst one among these three models in the 
viewpoint of ^ jdof . This is very different from the results of M3 obtained in Ref. [lB| where SCP 307 is used alone. 
For ACDM model, by using the datastes of SN la, BAG and GHD, we found the best fit values: f^^o = 0-283to °3i 
and /lo = 0.7321q;q44. And the corresponding minimum value of is Xmin = 327.246 {x^ /dof = 1.039). But, it is 
clear all the models give the same prediction that the universe is undergoing accelerated expansion and experience an 
transition from decelerated expansion to accelerated expansion. The transition redshift zt depends on the concrete 
parameterized form one assumed. Here Mi and M2 give the value about 0.6. However, the luminosity expansion 
model M3 give large transition redshift. From the right panel of Fig. [U one can not take the expansion model seriously 
for its bad behavior at relative high redshift, say z > 1. 

The 1(7 and 2(7 contour plots in qo — jo and qo — ho planes are plotted in Fig. [2] and Fig. [3l In qo ~ jo contour plots 
[21 the ACDM model result is included as a short line segment denoting la interval where the dot denotes the best fit 
value. It can be seen from Fig. [21 the best fit ACDM model is out of the range of Icr region in Mi and 2a region in 
M3, but it is in the range of la region in M2. It means that the best fit ACDM model is c omp atible at la and 2a 
confidence levels with M2 and Mi respectively. The results are compatible with that of Ref. [l6| except the model of 
expansion of luminosity distance. The results imply a universe with jo < ^1 from Mi and M2. 

IV. CONCLUSIONS 

In this paper, kinematic models are constrained by recent cosmic observations which include SN la, BAG and 
observational Hubble data (GHD). Here, we give three simple examples: Taylor expansion of q{a) at present oq = 1 
(Ml), constant jerk parameter j = jo (M2) and expansion of luminosity distance (M3). All of the models give the 
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0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 



z z z 

FIG. 1: The evolutions of deceleration parameters q{z) with respect to redshift z in la error regions. The left, center and right 
panels correspond to Mi {q{z) — qo + qiz/{l + z)), M2 (j — jo) and M3(expansion of luminosity distance) respectively. 




-1.5 -1.0 -0.5 0.0 -1.5 -1.0 -0.5 0.0 0.5 



FIG. 2: The contour plots of qo — jo with la and 2a regions. The left, center and right panels correspond to Mi {q{z) = 
qo + qiz/ (1 + z)), M2 {j = Jo) and Afs (expansion of luminosity distance) respectively, where the dots denote the best fit values 
of the parameters. The ACDM model result is included as a short line segment denoting la interval. 



same prediction of the evolution of the universe which is undergoing accelerated expansion at current and experiences 
a transition from decelerated expansion to accelerated expansion. But the transition redshift depends on the concrete 
parameterized forms of the models. The best fit values of the parameter of Afi and M2 predict the transition redshift 
is about zt ^ 0.6. However, M3 predict a large transition redshift Zt ^ 1. The /dof imply the same goodness of the 
models. But, from evolution curves of deceleration parameter as plotted in Fig. [l]and the knowledges of the history 



1.2 - 




-1.5 -1.0 -0.5 0.0 
9o 




FIG. 3: The contour plots of go — ho with la and 2a regions. The left, center and right panels correspond to M\ {q{z) = 
qo + qiz/ (1 + z)), M2 (j = Jo) and M3 (expansion of luminosity distance) respectively, where the dots denote the best fit values 
of the parameters. 
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of the universe, one found that Mi and M2 are better than M3 for the badness of behavior at high redshift of M3, 
say z > 1. 

The 1(7 and 2a confidence contours in qo — jo and qo — ho planes are plotted in Fig. [2] and Fig. [3] where the best 
fit value of ACDM model is denoted by a line segment, la interval. From the Fig. [21 one can find that Mi and M2 
are compatible with ACDM models at 2a and la confidence levels respectively. M3 is not compatible with ACDM 
model at 2a confidence level. For its badness of the evolution history of model M3 at high redshift. One can not take 
it seriously and treat it as a unreliable model. Then, one can conclude that the cosmic data favors the model having 
the value jo < —1 for the large parts of the confidence regions is under the line jo = —1- It is compatible with the 
conclusion obtained in Ref. [16] . 
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APPENDIX A: MODELS AND USEFUL RELATIONS 



In this appendix, one can find the analytical expressions between the Hubble parameter H{z), deceleration param- 
eter q{z) and jerk parameter j{z). 
Ml 



l + z 



(Al) 



H{z) = Ho{l + z)i+*+9i cxp ( - 



qiz 

l + z 



(A2) 



M2 



where 



J{z) = -qo - — l + z -77- — 72 ^ -r— 

l + z V (1 + z) l + z 



- 2 Uo 



l + z 



q{z) 



H{Z) = Ho[ci{l + ZY' + C2(l + z)"^]3 

^ Ci(l + Z)"^(^ - 1) + C2(l + zY-[^ - 1) 
Cl(l + Z)"i +C2(1 + Z)"2 

3{z) = jo 

1 

C2 a2 — 2 



Zt 



ci ai — 2 



1 
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(A3) 

(A4) 
(A5) 
(A6) 

(A7) 
(A8) 



2(l + go) -a2 , ^ 

ci = and C2 = 1 — ci 

ai — a2 



(A9) 



From (|A8p we see that jo < |. 

M3 - defined by the expanded luminosity distance Eq. ([5]), dL{z) = (z + Az'^ + Bz^), where A = (1 — qo)/'2 
and B = -{1 - qo ~ 3gg - jo)/6- 



H{z) = Ho 



[l + zy 



l + 2Az + {A + 3B)z^ + 2Bz^ 



(AlO) 
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l-2A~2{A + 3B)z- {A + 9B)z^ -2Bz^ 

'^^^^ " l + 2Az + {A + 3B)z^ + 2Bz^ ^ ' 

Jiz)^-[q + 2q^ + il + z)q'] (Al2) 

zt : the real root of 1 - 2A - 2{A + 3B)zt -{A + 9B)zf - 2Bzf (A13) 

ACDM, + r^A = 1 

H{z) = Ho + z)3 + (1 - rtm)] ^ (A14) 

qiz) = [(1 + zf - 2(l/r!„ - 1)] / [2(1 + zf + 2{l/n„, ~ 1)] (A15) 

j{z) = -1 (A16) 

zi = [2(f]-i-l)]^-l (A17) 
Note that the expressions for ACDM can be easily obtained from (I A4II A7|) . putting jo = —1 into 
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